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The unbinding properties of an interface near structured wedges are investigated by discrete mod- 
els with short range interactions. The calculations demonstrate that interface unbinding take place 
in two stages: i) a continuous filling-like transition in the pure wedge-like parts of the structure; 
ii) a conclusive discontinuous unbinding. In 2D an exact transfer matrix approach allows to extract 
the whole interface phase diagram and the precise mechanism at the basis of the phenomenon. The 
Metropolis Monte Carlo simulations performed in 3D reveal an analogous behavior. The emerging 
scenario allows to shed new light onto the problem of wetting of geometrically rough walls. 



I. INTRODUCTION 

Wetting phenomena concerns the properties of the 
liquid film which forms when an undersaturated vapor 
(fluid) is put in contact with a solid inert substrate. The 
critical properties of the liquid-vapor interface are, gen- 
erally, determined by nature and range of the intervening 
interactions. In this respect, e.g. for a planar substrate, 
different interactions potentials determine the applicabil- 
ity of partial or complete wetting regimes as well as the 
nature of the wetting transition (for a review on these 
phenomena see 0,0) 

However, substrate surface geometry can strongly in- 
fluence these properties. For example, in the complete 
wetting regime, adsorption isotherms of random rough 
substrates y, Q and linearly sculpted substrates [UIEQ] 
exhibit unusual exponents which are determined by the 
surface geometry. On the other hand, while for planar 
geometry partial wetting prevents the growth of macro- 
scopic films, in the same conditions, in pure wedges, we 
can have continuous filling phenomena by which the film 
thickness is driven to infinity [5|, |8( • Finally, and this is 
particularly relevant for wetting critical properties, there 
are many indications that increasing surface roughness 
can change the order of the wetting transition to first- 
order nEsmiii- 

In this paper we report the results of an accurate in- 
vestigation about the unbinding properties of a thermally 
fluctuating interface in wedge-modified systems (see Fig. 
1), we denoted structured wedges (SW). The choice of 
such a system structure is motivated by the fact that 
they incorporate contrasting geometrical motives (like 
wedges and ridges) which are typical of the geometry of 
rough surfaces. In this respect, while the separate effects 
of these geometries on wetting phenomena are already 
rather well known 0,^1' the effects of their combination 
were never approached in detail before. Here, we consider 
the study of both 2D (Fig. 1(a)) and 3D (Fig. 1(6)) 
SWs, and our results show that such structures implies 
two stage interface unbinding transition, one of which 
discontinuous. Our analysis of the phenomenon allows 
to understand its intimate connection with the surface 
geometry as well as to make some extension to the wet- 
ting properties of rough boundaries. 



II. THE MODEL 

In the framework of solid on solid (SOS) approach 
our model liquid-vapor interface corresponds to a lat- 
tice random walk (in 2D) or to a random surface (in 
3D) in the vicinity of a fixed substrate boundary. If 
Hx denotes the substrate boundary (integer) height at 
the position X (depending on 2D or 3D space dimen- 
sion, X denotes a single variable x or a couple of variable 
{x,y), respectively), the interface configurations can be 
specified in terms of the local relative (integer) height 
variables zx (see Fig. 2). The discrete nature of the 
model, implies that non horizontal walls (see the mag- 
nifying glass in Fig. 1) are shaped as staircases with 
slopes not exceeding unity. Consequently, the height 
function Hx satisfy the conditions H x+ i — H x = 0, ±1 
and H x+ i_ y — H xy = 0, ±1 in 2D and 3D, respectively. 
In 3D, because of the translation invariance along y axis, 
one has also H x ^ y+ i — H x y = 0. 




FIG. 1. Sketch of the 2D (panel (a)) and 3D (panel (&)) SW 
geometries studied in this paper. The magnifying glass shows 
the staircase nature of the tilted smooth walls (with slope 
|l/n|) by which the SW are constructed. 

At coexistence (no chemical potential differences with 
respect to its critical bulk value), an interface (see Fig. 
2(a), (6)), can be studied in terms of an Hamiltonian of 
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the form 
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where the first sum is done over all pairs of nearest neigh- 
bor columns. In the above expression £ and — U (with 
£,U > 0) are the energy cost of any interface step (or pla- 
quette in 3-D) and the energy gain of each interface con- 
tact (of horizontal step or plaquette) with the substrate, 
respectively; 7 = 1 or 00 determines the SOS or restricted 
SOS (RSOS) character of the implemented walk model, 
respectively. To limit the computational complexity, in 
3D we consider only RSOS interfaces. 
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FIG. 2. (a) Example of a configuration of a directed walk 
(heavy continuous line) with the ends anchored to a 2D 
boundary (light continuous line delimiting the shaded area). 
(b) An RSOS interface-surface configuration near our 3D SW. 



III. 2D STRUCTURED WEDGES 

In 2D, an interface in the vicinity of a boundary (with 
a generic shape) can be fully treated by a transfer matrix 
approach based on the matrices R x and L x defined as 
follows 
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Here u = e~ £ / kBT = e" 1 /* (* = k B T/£) and k = 
e u/k B T _ e «/i ( u _ \J j£) correspond to step and wall 
fugacities and, by definition, uj[R x ] z ^ z ' and w[L x ] ZyZ < are 
the Boltzmann weights of the elementary walks from 
(x, H x + z) to (x + 1, H x+1 + z') or (x - 1,14-1 + z 1 ), 
respectively. Thus, the partition function of, e.g. the 
walks between (xi,H Xl ) and (x2,H X2 ) (see Fig. 2(a)), is 
given by 
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But the transfer matrix approach allows also the di- 
rect calculation of the distance probability distribution 
function (PDF) at any position along the boundary. A 
walk like the one in Fig. 2(a) can be divided in two 
parts, e.g. at xq, obtaining two walks with Boltzmann 
weights [R Xl ^ Xo ] , z and [L Xo ^ X2 ] . z , respectively. Thus, 
the walk distance PDF at xq is proportional to the prod- 
uct [Rxi-tx ]o,z ■ [L Xo <~ X2 ]o, z . On the other hand, storing 
the initial walk distance distributions, at x\ and X2, in 
the vectors r Xl and l X2 (in the case of Fig. 2(a) we should 
set [r Xl ] z = [l X2 ]z = $z,o), the following iterations 



r x+ i = r x R x ; 



(4) 



allows to get the quantities [R Xl - 
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and [L a 



in terms of the components of r Xo and l Xo vectors. 



In- 
deed, iterations transfer toward right and left the lo- 
cal walk distance distributions to any position along the 
boundary, like in a forward diffusion process. Therefore, 
in the limit of an indefinite walk (our interface model), 
i.e. \xq — 2:1,2! ~~ > °°i the transferred right and left PDFs, 
defined as 
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will reach steady profiles depending only on the specific 
boundary conformation. In these conditions the interface 
PDF at any position xq will be given by 



and the quantity 
Af X0 = -\nP X0 (z) = 



(6) 



lnP^(z)+bxP^(z) , (7) 



can be seen as the corresponding local excess interface 
free energy profile. 

Our SWs boundaries contain only unit vertical steps. 
Thus, the transfer matrices can be of only three dif- 
ferent forms, i.e. those corresponding to a down step, an 
up step and no vertical step; we denote these matrices as 
D, U and F, respectively. For tilted wall regions with 
slope |1/ti|, the calculation of right and left PDFs can be 
done by implementing iterations Q in terms of F 71 " 1 !) 
or F n ~ 1 TJ matrices (depending on the descending or as- 
cending character of the wall); for large system sizes, the 
corresponding p( fl ) and p( L ) will asymptotically become 
eigenvectors of these matrices. In particular, in binding 
conditions (i.e. at high enough u values), the interface 
PDFs at x = ±w and x = (see Fig. 1(a)) will corre- 
spond to <l>^ and > 3\i and ty n being the bound eigen- 
vectors of matrices .F ra_1 .D and F n ~ 1 U, respectively. In 
this respect, an accurate investigation allows the deter- 
mination of these states (exact for n = 1,2, numerical for 
n > 2) both for SOS as well as RSOS walk models. 

A brief analysis of n = 1 case (tilted walls with unit 
slope) is particularly useful for the understanding of the 
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general scenario. Matrices F n ~ 1 D and F n ~ x U reduce to 
D and U and for SOS walks both and *S?i components 
scale as /if) and nfj , respectively. fj,£> 1 — , and 

thus the local character of $1 (i.e. \xn x < 1) is achieved 
until k > k{) i ci with k^ x c — u (^_^ ■ The correspond- 
ing interface PDF at the SW bottoms (i.e. x = ±iu) is 
P± w (z) = fi'jj and the resulting interface average dis- 
tance, (z)± w = % z 'p±Jti Z ) is § iven b y 



X/±W (2-k + k0J 2 )(l+2kLU 2 + k 2 LU i -k 2 UJ 2 )' w 

As k — > hp c , (z)± w diverges with continuity as 

[k — kD lt c] 1 an d thus ko ltC is the critical fugacity as- 
sociated to the continuous interface unbinding from SW 
bottoms. This transition is the analog of the critical fill- 
ing transition in pure wedge geometry |3, ll 1| . We have 
verified that for n = 2, 3,.. .,10 one observe the same 
behavior (also for RSOS walk model) and in Fig. 3 we 
show the corresponding phase-line boundaries, consisting 
in the the plot of the quantity u c (t) = t In ko n , c - The pic- 
ture demonstrate that in the wedge geometry, in contrast 
to the case of flat substrates (see the bottom curves in 
Fig. 3), a bound interface always (also at low t) requires 
a finite attraction energy u. On the other hand, the local 
minima in the phase-line boundaries confirm the possi- 
bility of reentrance phenomena which were already pre- 
dicted for rough self-affine boundaries. More details 
about these aspects will be reported elsewhere [l5|. 
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FIG. 3. SOS (continuous lines) and RSOS (dashed lines) 
interface phase-line boundaries in the u-t plane for 2D 
SWs (like in Fig. 1(a)) with tilted wall slope \l/n\ (n = 
1,2, ... , 10). Upper curves (n = 1, 2) are the result of exact 
calculations; the other curves were obtained by a numerical 
implementation of iterations (J1J. Lower (light) curves are the 
interface phase-line boundary for a flat wall [l6j|. 

On the contrary, fiu 1 = and therefore the lo- 

cal character of is achieved for k > kjj lyC with 



kjji.c = 1/(1 — uj) < &£>!,(;• Thus, when k approaches 
koi,c from above, and the interface unbinds from SW 
bottoms, it remains tightly bound to the central ridge. 
To be more precise, in the limit k — > [fc_D liC ] + the mean 
interface distance from the SW central ridge is finite and 
is given by 

M = ^(1-^ + ^) (Q) 

K /0 {I + U + U 2 ){1 + L0 2 ){\ -LU) 2 ' [ ' 

Finally, for k < kD llC we are out of binding condi- 
tions for D matrix and the lateral descending walls of the 
SW are no more able to bind the walks: iterations J3J, 
asymptotically produce delocalized right and left PDFs. 
In these conditions, the presence of the central ridge can 
have only very marginal effects (vanishing with the size 
of the structure) and, thus, the interface unbinds from 
the SW. The conclusion is: fco^c is also the threshold of 
a discontinuous unbinding transition! The curves in Fig. 
3 can be seen also as the interface phase-lines for 2D SW 
first-order unbinding. 
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FIG. 4. Excess free energy profiles (in arbitrary units) at 
x — ±w (panel (a)) and x — (panel (6)) for a 2D SW with 
wall slope \l/2\,d= 1000 and di = 2000 (see Fig. 1(a)). The 
free energy profiles have been obtained by a numerical imple- 
mentation of iterations Q at u = 0.2 (t = 0.6213 . . .). The 
curves are shifted vertically to increase the picture clearness; 
in panel (b) also a little horizontal shift has been adopted to 
avoid curves superposition at z = 0. 

The conclusive validation of the above, surprising, con- 
clusion is obtained by the analysis of the free energy 
profiles. In Fig. 4 we show a plot of Af± w (z) (panel 
(a)) and A/o(z) (panel (&)) profiles obtained at fixed 
t for an RSOS interface in a finite size SW with wall 
slope 1/2 (see the specific SW dimensions in the figure 
legend) characterized by an unbinding critical fugacity 
ko 2 ,c — 2.607536 . . . (exact calculation). Curves in panel 
(a) show, as k decreases, the following free energy pro- 
file evolution: i) a unique free energy minimum at z = 
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(for k — 2.6100 > ko 2 . c ) which delocalizes into a wide 

one at k = 2.6075 ~ fcl? SOS ' ) (he. the continuous transi- 
ts ,C V 

tion in the pure wedges); ii) double minima profiles for 
k = 2.6050,2.6025 < krj 2 , c with the minima placed at 
z = d± — d and z = d\ (i.e. coexistence between the 
state localized at the SW central ridge height and the 
bulk unbound state); Hi) dominance of bulk unbound 
state (fc = 2.6000 < fcc 3)C )- An analogous evolution is 
extracted by the free energy profiles of A/ in panel (6) . 
Therefore, in passing between ii) and Hi) we have a dis- 
continuous transition! Finally, the picture clarifies that 
the discontinuous nature of such transition is strictly re- 
lated to the fact that, for diverging system size, d pa- 
rameter maintains larger than the widths of the two free 
energy minima (scaling as ~ IT 1 / 2 ). 
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FIG. 5. Average interface distances (open dots) from 3D 
SW bottoms (panel (a)) and central ridge (panel (6)) and 
the corresponding distributions (gray shaded area). The data 
come from Metropolis Monte Carlo simulations at constant 
t — 2.0 and system sizes L = 100, W = 50, wall unit slope 
and di/d — 2 (i.e. d — 16). In the insets is the detail of 
interface distance distributions at the u values marked by the 
dashed lines. 



IV. 3D STRUCTURED WEDGES 

Transfer matrix approach is not applicable to the wet- 
ting problem in 3D. In this case we have considered 
Metropolis Monte Carlo simulations of an RSOS fluctu- 
ating interface (see Fig. 2(6)) in a 3D SW geometry like 
the one sketched in Fig. 1(6). Because of rapid increase of 
computation times with system size, we have considered 
relatively small systems: the simulations were done for 
squared system with size L — 50, 100, 150 (in the y axis 
direction) and W = L/2 (L and W are in lattice units) 
and using periodical boundary conditions along both x 



and y axes. The SWs considered in the simulations were 
usually constructed by (staircase) tilted walls with slope 
1 or 1/2 and d\/d ~ 2 (see Fig. 1(6)). As usual, the im- 
plemented Metropolis Monte Carlo algorithm was based 
on local moves corresponding to the attempt of chang- 
ing (at random) the local surface height by one unit (i.e. 
z x .y - > z x ,y ± 1). After equilibration, the calculation 
of equilibrium average parameters has been performed 
on the basis of very long simulations of up to 10 5 -j- 10 6 
MCS (1 MCS ee L 2 Monte Carlo moves). 

In Fig. 5, as a representative of the general behav- 
ior, we show the u dependence of the interface averages 
distances from SW bottoms ((z)± w , panel (a)) and SW 
central ridge ((z)o, panel (6)) obtained a constant t for 
a given 3D SW (see geometrical specifications in figure 
legend). The behavior of these two quantities follows a 
scheme very similar to the one outlined for 2D SW in- 
terface unbinding: i) a continuous filling-like transition 
of the two component wedges at u ~ 1.607 (note also in 
the right inset in panel (a) the roughly flat height dis- 
tribution); ii) a discontinuous detachment from the SW 
central ridge at u ~ 1.309 (see the double peak struc- 
ture of the interface height distributions in the left insets 
in panel (a) and (6)). However, at least for the finite 
size system considered here, continuous and discontinu- 
ous unbinding seem to be separated by a finite u gap. 
Most probably this is due to the linear extension of the 
system along y axis which corresponds to the main differ- 
ence between the present 3D SW geometry and the 2D 
one. With an extensive simulations program, we planned 
to evaluate the whole interface phase-lines for 3D wedges 
with the intent to clarify also this aspect. 



V. CONCLUSIONS 

Our investigation demonstrate that interface unbind- 
ing in SWs takes place in two stages: the first one corre- 
sponds to the continuous filling transition of the two com- 
ponent wedges; in the second stage the interface jumps, 
discontinuously, from a state localized at the SW cen- 
tral ridge and the bulk unbound state. We stress the 
fact that in SWs the bulk unbound state come into field, 
rather than in some artificial way p"?|. just because of the 
corrugated geometry of the substrate. First-order nature 
of the final unbinding transition is due to the competi- 
tion between binding at the central ridge and whole SW 
filling. In other words, it is the spatial combination of 
wedges and ridges of SWs to create the conditions for a 
discontinuous interface unbinding. As seen in section III, 
a single ridge is only able to strongly pin the interface 
retarding (with respect to a wedge) its unbinding from 
the apex (for more details see also ref. [l5|h In complete 
wetting regime, as shown in ref. [T^ | , interface unbinding 
near the apex mimics planar critical wetting. 

Generalization of our conclusions to more complex 
SWs and/or to random rough surfaces should take into 
account, carefully, the type and the scaling properties of 



■5 



the corrugation introduced by the specific boundary. On 
the other hand, our analysis about interface unbinding 
in 2D SWs give some more efforts on the key role of sur- 
face roughening exponent Cs in determining the nature 
of the wetting transition in self-affine rough substrates. 
In the light of the present results, we have started new 
investigations of 2D and 3D rough geometries, with the 



purpose to get conclusive insights in this issues. 
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